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The effect of a ring exchange on doped two-leg ladders is investigated combining exact diagonal- 
ization (ED) and density matrix renormalization group (DMRG) computations. We focus on the 
nature and weights of the low energy magnetic excitations and on superconducting pairing. The 
stability with respect to this cyclic term of a remarkable resonant mode originating from a hole pair 
- magnon bound state is examined. We also find that, near the zero-doping critical point separating 
rung-singlet and dimerized phases, doping reopens a spin gap. 

PACS numbers: 74.20.Mn, 75.10.-b, 75.40.Mg 
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Two-leg cuprate ladders 1 and two-dimensional (2D) 
high-temperature (high-T c ) superconductors seem to 
surprisingly share many similarities. In addition, ladders 
exhibit exotic physical behaviors characteristic of one- 
dimensional systems. Under high hydrostatic pressure, 
the intrinsically doped Sr2Cai2Cu2404i+,5 ladder mate- 
rial 2 becomes superconducting hence reinforcing the sim- 
ilarity between two-leg ladder and layer-based cuprates. 
Doped two-leg ladders are also appealing for theorists 1 : 
its ground state (GS) is believed to be a simple realization 
of the Resonating Valence Bond (RVB) state proposed 
by Anderson 3 in the context of high-T c . In addition, 
ladders are one of the simplest realistic correlated elec- 
tron systems to investigate novel mechanisms of super- 
conductivity 4 . Interestingly, new low-energy magnetic 
excitations have been observed by Nuclear Magnetic Res- 
onance (NMR) under pressure, i.e under doping the lad- 
der planes 2 , a finding possibly connected to predicted 
exotic magnon-hole pair states 5 . 

Taking into account ring exchange has proven to be 
essential to understand the spin dynamics of insulating 
cuprate materials such as both ladders 6 and high-T c 7 . 
Ring exchange which permutes cyclically spins on square 
plaquettes is known to suppress magnetic order in sev- 
eral 2D systems. It naturally appears in a fourth order 
expansion in t/U of the Hubbard model near half- filling 
and also in ab initio calculations 8 . In order to study the 
properties of lightly doped ladders 9 with ring exchange, 
we shall use the framework of the t — J Hamiltonian and 
include a cyclic exchange of magnitude K 
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with Vg Gutzwiller projectors. Labels i,j,k,l of per- 



mutation operators Pn^ stand for the sites encountered 
along each square plaquette. The expression of the cyclic 
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includes bilinear frustrating terms and biquadratic 
terms 10 . The K term leads to a rich phase diagram at 
half- filling 11 . In particular, a phase transition between 
the RVB rung-singlet phase and a phase with staggered 
dimer order occurs 10-13 at the critical point K c /J ~ 0.2. 
Experimentally proposed 6 values of K / J for the insulator 
ladder material range in a window 0.025 < K/J < 0.075. 
We are also interested here in the proximity of this crit- 
ical point. Thus, we shall restrict ourselves to param- 
eters —0.2 < K/J < 0.2 and assume a physical value 
of J/t — 0.5 for which superconducting fluctuations are 
dominant at low doping 9 when K/J = 0. A similar 
model has recently been studied on the 2D square lat- 
tice 14 in a large N limit of a derived SU(N) model. 

It is instructive to consider in more detail the effect of 
the cyclic exchange in terms of the RVB description of 
ladders 15 . Consider a single plaquette. With K = 0, the 
ground state is exactly described as an RVB state, with 
the horizontal dimer state added to the vertical dimer 
state 16 . A small nonzero K rotates between horizontal 
and vertical valence bond states and acts either to en- 
hance (K < 0) or suppress (K > 0) the resonance. The 
ground state has energy E = — 3J + 2K. For K = J/2, 
the resonance is destroyed and the singlet ground state 
is doubly degenerate between linear combinations of the 
two VB states. On the other hand, a nonzero K of either 
sign is frustrating to Neel order. Thus the cyclic term 
mediates between VB and RVB order. 

In this article, after briefly reviewing the undoped case 
in connection with previous work, we first consider the 
case of two holes in a ladder, a simple limit where we 
expect to get physical insight on the various elementary 
excitations carrying different spin and charge quantum 




FIG. 1: (Color online) (a,b) Spectra for undoped ladders 
from ED calculations on ladders of length L — 12, 14 and 
16. Grey areas are guide to the eyes, (c) Spectra for ladders 
doped with two holes and L = 9, 11, 13. Note that k x — n 
does not belong to the Brillouin zone of systems with odd L. 
The dashed line shows A p computed independently. 



numbers. In a second step, we address the finite doping 
issue where we identify the previous excitations, discuss 
the effect of doping and connections to experiments. 

Elementary excitations - The magnetic energy spec- 
trum is shown in FIG. 1 (a). One can distinguish a 
one magnon branch excitation around its minimum at 
momentum k = (n, tt). The bottom of a two-magnon 
continuum is located around k = (0, 0) with an energy 
about twice the spin gap. When K/J increases, the 
spin gap is reduced and vanishes at the critical point 



FIG. 2: (Color online) Comparison between characteristic en- 
ergies of elementary magnetic excitations in doped ladders. 
Data are extrapolated from DMRG computations on systems 
of length up to L — 48. Proposed experimental values of Kj J 
for undoped systems are in the Exp. window. 



as was previously proposed. There (FIG. 1 (b)), the 
triplet branch is divided into a symmetric mode (k y = 0, 
k x G [0, 7r/2]) reminiscent of the two-magnon continuum, 
and an antisymmetric mode (k y = tt, k x £ [tt/2, 7t]) rem- 
iniscent of the one-magnon branch. The doped system 
(FIG. 1 (c)) is a Luther-Emery liquid provided K/J re- 
mains small enough : a gapless charge mode appears at 
low momenta, corresponding to the motion of the hole 
pair, while the one magnon branch remains gapped. A 
quasiparticle (QP) continuum starts for an energy which 
equals the pairing energy (computed independently, see 
after) and larger than the spin gap. Note that the triplet 
branch has its minimum at an incommensurate momen- 
tum k d - = (7r(l — 5), tt) and that, in the case of two holes 
under study, k,5 — > (tt, tt) in the thermodynamic limit. 
As K/J increases, the spin gap is reduced as is the onset 
of the QP continuum. Lastly, the system with K < has 
similar behavior to a two-leg ladder with a rung coupling 



J± > J 



In order to distinguish between the collective magnon 
and the QP excitations, we first focus on ladders with 2 
holes forming a singlet d-wave bound state when L — > oo. 
Following ref. 6, we define 

A M = E(0h,S=l)-E(0h,S = 0) 
A s = E(2h, S = 1) - E(2h, S = 0) 
A p = 2E(lh,S = 1/2)- E(2h., S = 0)- E(0h,S = 0), 

where E(nh, S) is the GS energy with n holes and spin 
S. Am represents the magnon gap, A 5 the spin gap with 
two holes, and A p the pairing energy. Each gap shown 
on FIG. 2 has been calculated on finite systems with 
DMRG 17 and then extrapolated to the thermodynamic 
limit. First, let us make a few general statements on these 
energies : (i) since a magnon can be created arbitrarily 
"far away" from the holes, then A5 < Am, (ii) breaking 
the hole pair creates two quasiparticles so that A 5 < A p . 
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FIG. 3: Hole-hole static correlation functions from ED com- 
putations on L — 13 ladder with two holes for different K/J. 
The area of circles is proportional to the correlation value. 
Up left circle indicates the mean expectation value. 



f 0. 















r— 









K/J = -0.08 
w = 12.8% 

K/J = 
w = 1 1 .8 % 

K/J = 0.12 
w = 9.5 % 



FIG. 4: Spin density around the hole pair in the first triplet 
state from ED data on L — 13 ladder with two holes for differ- 
ent K/J. Weights w of the projected function are indicated. 
Lengths of arrows are proportional to the local density values. 



These inequalities are well verified in the data. The dis- 
continuity of the spin gap in the infinitesimal doping limit 
has been attributed to the existence of a bound state be- 
tween the hole pair and the magnon, responsible for a 
sharp resonating mode at finite density similar to what 
can be observed in Inelastic Neutron Scattering (INS) ex- 
periments on underdoped cuprates. The binding energy 
of this bound state can be defined as min(A p , Am) — As, 
see FIG. 2. When A p < Am, the continuum starts at 
the QP continuum so that the magnon is scattered by QP 
and acquires a finite life time (dashed part of the Am line, 
K/J < —0.04). When Am < A p , we expect a contin- 
uum of hole pair-magnon scattering states above Am and 
below the QP continuum, as shown by the grey dashed 
area for —0.04 < K/J < 0.2. The binding energy of the 
bound state is strongly affected by the cyclic exchange 
with an instability treshold around a value of K signifi- 
cantly smaller than the one of the critical point. More- 
over, the pairing energy vanishes simultaneously with the 
spin gap when K reaches K c . 

The binding of holes (and also spinons) in the two-leg 
ladder can be described as an effect of the resonance in 
the RVB language 15 : an isolated hole induces a stag- 
gered, non-resonating VB state either to the left or the 
right of the hole, which is only healed by another hole. 
The destruction of resonance as K is increased towards 
K/J w 0.2 should translate directly into weakened and 
eventual loss of pairing. This effect is clearly evident in 
FIG. 3. For K < 0, the pairing is strengthened both by 
the enhancement of resonance and by the fact that the 




FIG. 5: (Color online) Spin gap for doped two-leg ladders as 
a function of hole doping 8. Data extrapolated from DMRG 
computations on systems of length L = 16 to 64 and doping 
8 = n/16, n — 0, • • ■ , 6, with a number of kept states up to 
m = 1600. The 8 — magnon energies are indicated by grey 
symbols (see arrows). Half the pairing energy for a finite 
system is also given by dashed lines. 



cyclic exchange does not act on a plaquette with a hole, 
favoring states with holes sharing plaquettes. 

The pair-magnon binding can be viewed as closely re- 
lated to ferromagnetic spin polarons. The hopping of 
holes in a ferromagnetic background is unfrustrated; the 
magnon produces a small region of enhanced ferromag- 
netic correlations. In FIG. 4, the bound state has been 
projected onto the states with two holes on the diagonal 
of a plaquette (of weights w given on figure) showing in- 
deed a strong ferromagnetic region for K < 0. More pre- 
cisely, the hole pair-magnon binding results from a gain 
in kinetic energy coming from processes that exchange 
the triplet and the pair locally on a plaquette 18 . As K/J 
increases, the "free" magnon excitation lowers its energy 
while the above resonant processes are less affected. 

Finite doping study - We now turn to the more com- 
plex case of finite doping. As shown below, we find signa- 
tures of the various excitations identified in the previous 
part. The spin gap for several hole density 8 is displayed 
in FIG. 5. For K/J = 0, we see that it is reduced un- 
der doping because holes weaken spin exchange 19 but re- 
mains robust 20 . For small magnitude of cyclic exchange 
(for instance K/J = 0.06) this effect is still significant, 
leading to an appreciable lowering in comparison with 
the undoped system. As the critical point is approached, 
the magnetic parent acquires a smaller spin gap but the 
relative decrease under doping is less and less significant. 
Interestingly, at the half-filling critical point K c /J ~ 0.2, 
doping eventually leads to a reopening of the spin gap. 
Qualitatively, this could be explained by the fact that a 
hole on a plaquette prevents the action of the ring ex- 
change. Therefore, at sufficiently large doping (S ~ 0.4), 
the spin gap depends weakly on K for studied values. 
Pairing energy also increases at large densities which sug- 
gests a deep relation between pairing and spin gap. 

We have computed static and dynamic spin structure 
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FIG. 6: (Color online) (a) static structure factor <S(k) for a 
ladder of length 12 with 4 holes (8 = 0.17). (b) 5(kj)/ <<S) for 
different doping as a function of K/ J, from systems with L 
ranging from 8 to 16. 5(k^) is normalized w.r.t. the average 
structure factor (S) = (2L)" 1 £ k S(k) = |(1 - 5). 
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FIG. 7: (Color online) (a) dynamic structure factor <S(k,(jj), 
with k = (k x , k y ), for a ladder of length 10 with 4 holes. Red 
dots indicate poles and their relative weight (see text), (b) 
relative weight of the resonant state at k*. 



factors for undoped and doped systems by ED with stan- 
dard continued-fraction techniques and studied their evo- 
lution with increasing K/J. The static structure factor 
<S(k) gives insights on possible magnetic ordering. The 
dynamic spin structure factor <S(k, u>) probes magnetic 
excitations that can be detected in INS experiments. We 
use the standard definitions 

5(k) = (S*_ k Sl) = J dre^ (S*(0)S*(r)) . 

5(k,w) = --3(Vo|g k ,„]. ^h/>o), 

it uj + h,Q + ir] — n 



where Eg is the ground state energy and r\ gives a finite 
broadening to the peaks. Within our choice of normal- 
ization, iS(k^) should either diverge as L a , a > or lni 
for quasi-ordered states or converge for states with fast 
decaying correlations. At fixed L, S(ks) follows the fi- 
nite correlation length and shows a maximum defining 
a critical point K C (L). In the thermodynamic limit, 
K C (L) — > K c . Such a behaviour is clear at half-filling 
(see FIG. 6 (b), black circles corresponding to L = 12, 14 
and 16). From the data at finite density, the incommen- 
surate nature of the resonant mode is clearly seen and 
consistent with the expected value k^. This feature is 
preserved up to the critical point. Moreover, in contrast 
to the undoped case, the maximum of 5(k^) with K is 
hardly noticeable at large density. This is consistent with 
the finite value of the spin gap found at large doping. 

Inspecting <S(k, u>) reveals a pole well separated from 
the rest of the spectrum at momentum k = k^ (FIG. 7 
(a)). This peak can be clearly identified as the resonant 
mode issued from the triplet BS discussed above 5 . Al- 
though ED suffers from finite size effects and dynamic 
spectra are tedious to analyze due to multi-excitation 
continua, the evolution of its relative weight given by 
| (^Re S | |V>o) | 2 /«S(k 5 ) (where |V>Res) is the resonant 
state wave function) is still very instructive. Note that in 
the thermodynamic limit, the resonant mode should not 
give rise to a single pole but rather a singularity 20 . The 
evolution of the relative weight is displayed in FIG. 7 (b) 
and shows that the BS remains robust at finite density 
and K > 0, keeping a sizeable fraction of the total weight. 
Consequently, we predict that it should be observable 
in spin dynamics sensitive experiments. This conclusion 
on the robustness and observability of the Luther-Emery 
phase is in agreement with recent NMR experiments on 
doped ladder materials 21 . 

Conclusion - In summary, we show that the nature 
of the Luther-Emery phase and its remarkable resonant 
mode is preserved for intermediate physical K/ J values. 
In this regime, both spin gap and superconducting pair- 
ing energy are very sensitive to the ring exchange. In 
addition, doping is shown to have a drastic effect at the 
magnetic critical point of the parent Mott insulator. 
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